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g ; 1 Introduction 

Let be a probability distribution on a finite dimensional linear space E depending on 
^ ' a parameter A which belongs to a subset A of M^. Suppose that 

/"A = f{x,\)(^{dx), 

where a is some reference measure, and that for each x in E, the map A i— > /(x. A) defined 
on A is measurable. For a probability distribution u{dX) on the set A, define 



h{x) = / f{x,\)v{d\). 
J A 

Then the probability measure 

fiuidx) = h{x)(j{dx) 

is called the mixture of the distribution ii\ with respect to v. (See Feller [4], Vol. II, page 
53 or Johnson et al. [8], page 360). Usually, ^\ is said the mixed distribution and v the 
mixing distribution (see Karlis and Meligkotsidou [9]). 



"Corresponding author. E-mail address: abdelhamid.hassairi@fss.rnu.tn 



1 



A special case of interest is when is not concentrated on an affinc hyperplane of E, ^\ is 
the A— power of convolution of ji, and A is the so called Y0rgensen set of Specifically, 
let 

L^{e)= I cxp{{e,x))n{dx) (1.1) 

Je 

denote the Laplace transform of where (,) is the duality crocket, and suppose that the 
set 

= interior{0 G E*; 1^(9) < +00} (1.2) 

is nonempty. Then the set 

A = {A > 0; 3 nx such that L^^ {0) = {L^{9)f for aU 9 G G(/i)} (1.3) 

is called the Y0rgensen set of ^ and the measure ^\ is its A— power of convolution. This 
parameter A appears in the most common models, it is in particular, the intensity in a 
Poisson model, the variance in a Gaussian model, and the shape parameter in a gamma 
model. For A and A' in A, we have that 



Ma * /"A' = Ma+a' • 



The set A contains always the set IN* of positive integers and it is equal to ]0, +00 [ if and 
only if ^ is infinitely divisible (sec Scshadri [13], page 155). 

In fact, if Xi, . . . , Xj^ are iid random variables with distribution ji, then the distribution 
of Xi + . . . + Xjv is /xat = A**'^, the iV— power of convolution of ji. Accordingly, for 
any distribution ^ and any positive integer iV, one may consider the distribution as 
defined in (1.3). When ^ is discrete, i.e., with countable support, then the mixture of ^ 
with respect of a distribution u on the parameter A'' is known as a compound distribution. 
The most famous compound distribution is the one corresponding to the case where u is 
Poisson (see Feller [4], Vol. I, page 286 or Vol. II, page 451 or Aalen [1], or Y0rgensen 
[14], page 140). In fact, the real Poisson distribution appears in numerous works either as 
a mixed distribution (see Bhattacharya and Holla [2] or Johnson et al. [8], page 366) or 
as a mixing distribution (see Per line [12]). In the present work, we will be interested in a 
very special case in which the mixed distribution is defined on the cone of (r, r) positive 
definite symmetric matrices Q with a parameter which belongs to a subset of IR*^. More 
precisely, the mixed distribution will be the absolutely continuous Riesz model introduced 
in Hassairi and Lajmi [6] 
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This model contains the Wishart model, since i?(s, cj) reduces to a Wishart distribution 
when si = S2 = . . . = Sd and it has a convolution property which is analogous to the one 

satisfied by the ordinary powers of convolution. In fact, if s and s' are in — - — , +00 

i=i-l ^ 
then 

R{s, a) * R{s', a) = R{s + s', a). 

The mixing distribution will be the multivariate Poisson distribution on W. For simplicity, 
we will be interested in the case where a is equal to the identity matrix of size r denoted 
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Ir- We first show that the mixture distribution is expressed in terms of the modified 
Bessel function. We then determine the domain of the means of the generated natural 
exponential family, and we calculate its variance function. This provides a rich class of 
natural exponential families. 

2 The Riesz exponential dispersion model 

A dispersion model is a class of natural exponential families where each family is a power of 
convolution of the other. In this section, we will first recall some general facts concerning 
the exponential dispersion models in an Euclidean space. Then we introduce the Riesz 
dispersion model on symmetric matrices. 

2.1 Exponential dispersion model 

Let E be an Euclidean space with finite dimension n, and let (, ) denote the scalar product 
in E. If ^ is a positive measure on E, we denote by Ai{E) the set of measures fi such that 
0(//) given in (1.2) is not empty and /j, is not concentrated on an afline hyperplane of E. 
If /X is in M{E), we denote 

k^{e) = logL^{e), for all 9 in e(/x) 

the cumulant function of /x, where is the Laplace transform of /x defined in (1.1). 
To each /x in M{E) and 9 in 6(/x), we associate the probability distribution on E 

P{9,n){dx) = exp{{9,x) - kf,{9)) nidx). 

The set 

F = F{i,) = {P{9,fiy, 9ee{fi)} 

is called the natural exponential family (NEE) generated by /x. We also say that /x is a 
basis of F. 

The function is strictly convex and real analytic. Its first derivative A;^ defines a 

diffeomorphism between 0(/x) and its image Mp. Since A;|^(0) = J xP{9, fi){dx), Mp is 

called the domain of the means of F. The inverse function of k'^ is denoted by '0^ and 
setting P{m,F) = P('^^(m),/x) the probability of F with mean m, we have 

F = {P{m,F);me Mp} , 

which is the parametrization of F by the mean. 

Now the covariance operator of P{m,F) is denoted by ^^(m) and the map 

Mf LsiE); m ^ Vpim) = k'/^i^m)) 
is called the variance function of the NEE F. It is easy proved that for all m G Mp, 

Vp{m) = (V;(m))-\ 

and an important feature of Vp is that it characterizes F in the following sense: If F 
and F' are two NEFs such that Vp{m) and Vp'{m) coincide on a nonempty open set of 
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Mi?nMi?/, then F = F' . In particular, knowledge of the variance function gives knowledge 
of the natural exponential family. 

Let II be an element of Ai{E) and let A be its Y0rgensen set defined by (1.3). Then the 
set 

{P{e, A) = exp {<e,x> -kf,^ {9)) n; ee e{n), A G A} 
is called the dispersion model generated by ji. For more details, we refer to Letac [11]. 

2.2 Riesz natural exponential families 

Let E be the Euclidean space of (r, r) real symmetric matrices equipped with the scalar 
product {x, y) = tr{xy), and the inner product x.y = - {xy + yx^ , where xy is the ordinary 
product of two matrices. 

We denote by ei, 62, . . . , the canonical basis of If; = (0, . . . , 0, 1, . . . 0), (1 in the 

i^^ place), and we set q = diag{e.i) for all 1 <i <r. 

For X & E, we consider the endomorphism L(x) of E defined by 

L{x) : y\ — > x.y 

and we set 

P(x) = 2(L(x))' -L(x'). 

We denote by f2 the cone of (r, r) real symmetric positive definite matrices. 
For X = {xij)i<ij<r in E and 1 < A; < r, we define the sub-matrices 

Pkix) = {Xij)i<ij<k and Pfc(x) = {Xij)r-k+l<i,j<r- 

For convenience, Pk{x) and P^{x) are also considered as elements of the space E, where 
the other entries arc equal to zero and we set Pq{x) = 0. 

Let Afc(a;) and A^(x) denote the determinant of Pkix) and the determinant of P^ix) 
respectively. Then the generalized power of a; in J7 is defined, for s = (si, S2, • • • , Sr) £ , 

by 

A,(x) = Ai{xy'-'^A2{xy^-'^ ArixY^ (2.4) 

Note that if for all i G {1, . . . , r}, Sj = p G M, then As{x) = (deta;)^'. We also define 

Atix) = {Alix)r-^-{A*,ix)r-^^ . . . (A;_l(x))^^--^^(A;(x))^^ (2.5) 

It is shown (see Hassairi and Lajmi [6]) that for all x G and all s G M'', we have 

A,(x-i) = Al,,(x), (2.6) 

where s* = (s^, Sr-i, ■ ■ ■ , si). . 

We denote by Tj^ the set of lower triangular matrices with positive diagonal elements, 
and for u G Tj^, we define on E the automorphism 

u{y) = uyu*, 

where u* denotes the transpose matrix of u. 
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It is well known that for all x G fi, there exists a unique u G T^^ such that 

X = u{Ir), 

where is the identity matrix of order r, it is the Cholesky decomposition of x. 
We also have (see Hassairi and Lajmi [6]) that for all 1 < z < r, 



and for all s = {si,S2, ■ ■ , Sr) G JEC, 

As{u{Ir)) = A*_Au*-Hlr)). 

Recall also that for x G and u G T^^, we have 

Ai{u{x)) = Ai{uiIr))Aiix) = ul... u^Ai{x), 
where for alH G {1, . . . , r}, 

Ui = {u,Ci). 

(See Faraut and Koranyi [3], page 114). 

Now let E be the set of elements s = (si, S2, . . . , s^) G defined as follows: 
Consider the function ^ defined from IR"'' into {0; 1}, by 



For all {ui,U2, ■ ■ ■ , Ur) G ]R!j_, we define 

Sl = Ul 



ifa = 0, 

1 if a > 0. 



(2.7) 

(2.8) 

(2.9) 
(2.10) 



,^^ e(^i) + - + e(^.-i) ^VfcG{2,...,r}. 



A result due to Gindikin [5] and proved in Faraut and Koranyi [3], page 124, says that 
there exists a positive measure Rg such that for all 9 G — ri, 



LRXe)= I e(''-^Rs{dx) = As{-e-') 
Je 



if and only if s is in H. This measure is called the Riesz measure with parameter s. 

i-l 



When s = (si, S2, . . . , Sr) € S \ 



i=l 



-, +00 



, the measure Rg is concentrated on the 



boundary dVi of Vt and when s = (si, S2, . . . , s^) is such that for all i, Sj > 



i - 1 



the 



measure Rg is absolutely continuous with respect to the Lebesgue measure and is given 
by 

Rs{dx) = ^ K{x){dx) 
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r(r + 1) 

with n = ^ — - the dimension of E and 



Tnis) = {27r)^ Y[t(s, - ^ ■ (2.11) 



i=l 



It is shown in Hassairi and Lajmi [6] that the measure Rg generates a natural exponential 
family if and only if si 7^ 0. In this case 

F = F{Rs) = i^R{s,a) = ^^-^R,, a e s € S, si 7^ | . 

The distribution R{s,a) is called the Riesz distribution with parameters s and a. It is 
shown in Hassairi and Lajmi [6] that the Laplace transform of R{s,a) is defined for 6 in 
o" — ri, by 

LR{s,a)[B)- — ^^^^,1^ — • (2.12) 
This implies that if a is an element of and if s and s' are in S, then we have 

R{s, a) * R{s', a) = R{s + s', a). 

i — 1 

Let a be an element of Q. If s satisfies the conditions Sj > — - — , for 1 < i < r, then the 
Riesz distribution is given by 

e-<^'^>A,_n(x) 

When si = S2 = ■ . . = Sr = p > 0, R{s, a) reduces to the Wishart distribution with 

r — 1 

parameters p > — - — and a G fi, 

Wip,a)idx) = r^(^) det(a-P) '"'^'^'^^ dct(x)^--7 1, (x) (dx), 
with Laplace transform equal for all ^ G cr — fi, to 

3 The mixture of the Riesz distribution with re- 
spect to the multivariate Poisson 

Consider the Poisson distribution on IN^ with parameter A = (Ai, A2, . . . , A^) G (M"'")'' 

V^r A^ 

u{dx) = e-2^i=i^' — 6q{dx), (3.13) 
where q\ = qil q2\ ■ ■ ■ qr^. and A'^ = A^^ ^^'2 ■ ■ ■ ■ Then for all 6 G IR*^, 

r 

L^e) = JJe^'("''-i). (3.14) 

i=l 
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The variance function of the natural exponential family generated by u evaluated for 
m g]0, +oof is equal to 

V{m) = diag{m). 

Note that this variance function is homogeneous of order 1 and consequently the natural 
exponential family F{v) belongs to the Tweedie scale of M'' (see Hassairi and Louati [7]). 

Next, we will show that the mixture of the distribution Ir) with respect to v is related 
to the modified Bessel function of the first kind and of order 1. 
Consider the second-order linear differential equation 

z''y"{z) + zy'{z) + {z"" - a^)y{z) = 0, (3.15) 

where ^ is a complex variable and a is a real parameter. One of the solution of (3.15) is 
the function J{a, z) known as the Bessel function of the first kind of order a, and defined 
for 2; in (C by the series 



(-i)'(l) 



2k+a 



In many applications, one frequently encounters two functions /(a, z) and K(a, z) called 
respectively the modified Bessel functions of the first and second kinds of order a, defined 
on the complex plane cut along the negative real axis by 

I{a,z) =r''J{a,iz), (3.17) 

and 

n I(— a, x) — /(a, x) ,„ 

K{a,z) = -^ ; ' ^ a^0,±l,±2,.... 

2 sm(all) 

It is an immediate fact that I(a, z) and K(a, z) are linearly independent solutions of the 
differential equation: 

z'^y"{z) + zy\z) - {z^ + c?)y{z) = 0. (3.18) 

For more details about the Bessel functions, we refer to Lebedev [10], page 108 and to 
Feller [4], Vol. II, page 58. 

We also need to establish the following result related to these functions. 
Proposition 3.1 Let x > and let b> 0, then 



Proof 

Consider the function 



Deriving gives 
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This impUes that 

fcelN 

Consequently, gf,{x) is a solution of the differential equation 

xY{x) + hx^-^y\x) - x^-\{x) = 0. 
Or equivalently, g^^ {x) is solution of the equation 

xy"{x) + hy'{x) = y{x). (3.20) 

Using now the fact that /(a, z) and K[a, z) arc two linearly independent solutions of (3.18), 
we deduce that x~ I{b — 1, '2^/x) and x~ K{b — 1, 2^/x) are two linearly independent 
solutions of (3.20). It follows that the general solution of (3.20) is of the form 

cix^ I{b - 1, 2Vx) + C2X^ K{b - 1, 2Vx), 

where ci and C2 are two constants. 

According to (3.19), (3.16) and (3.17), we deduce that 

Therefore ci = 1 and C2 = 0. 
Consequently, for all x > 0, 

□ 



3.1 The probability density of the mixture of the Riesz dis- 
tribution with respect to the multivariate Poisson 

Let p= ^0, ^, . . . , ^-f^^. A; = (fci, k2,-- - , kr) be in W and let 

Rk = R{k + p,Ir)- 

Suppose that k = (ki,k2, ■ ■ ■ ,kr) has the multivariate Poisson distribution ly defined in 
(3.13). For simplicity, in what follows, we will denote by p the mixture of Rk by u. The 
following theorem gives the expression of p in terms of the modified Bessel function of the 
first kind. 

Theorem 3.2 

^^idx) = X-Tl ri " f\ I (l,2yA,Ae,(x)) ljx)(dx), 

(27r)^V^^»=i VAi-i(x) V " J 

where Ao(a;) = 1. 
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Proof 

Let a = (ai, 02, . . . , Or) be an element of ]0, +oof and let 

Rk,a = R{k + p + a,Ir). 
Denote the mixture of Rk,a with respect to v. Then 

IJ,a{dx) = ha{x)l^{x){dx), 

where 



ha{x) = e-^^=.'^ Y: ,^ , 2 \ ' (3.22) 



and 



ge IN 

According to (2.4) and (2.11), (3.22) may be written as 

(2n)^ Vdit(i) »ii V 'i'- ^(^^ + V Ai_i(x 



Therefore 

ha{x) = n ^rT=^ ^e.(^)"'5a. (A,Ae,(x)) , (3.24) 

(2n) 4 v^d^t(^i=i VvAi(a;) y 

where ga^ is defined in (3.19). Inserting now (3.21) in (3.24), we obtain 

ha{x) = n ''/A / ^ Ae,(x)^/ ai-l,2^AiAe,(a 

(2n) ^ ydit(i) »=i \ V A» (x) V V 

According to (3.23), we deduce that 

p-tr{x) r ( /TT „-Aj / ^ \ \ 

= ^ n / -l,2^A,Ae,(a;) K{x){dx). 

(2n)^ydit(i)«=i VvA^-iW V V yy 

Invoking the fact that for all x > 0, we have 

/(l,x) =/(-l,x), 

(see Lebedev [10], page 110), the proof of the theorem is complete. □ 
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3.2 The variance function of F{ii) 

In this subsection, we study the natural exponential family F generated by /x. We first 
give the Laplace transform of n, then we determine the domain of the means and the 
variance function of the family F. 

r 

In what follows, for r > 1, we denote = ^ - e^- and we set = 0. 

Theorem 3.3 For all 6 e Ir — we have 

L,{e) = Al_^{Ir - 9) exp (j2 Xi (a1,^_,^,(/. - ^) - l))- (3.25) 



u=i 



Proof 

Let Xk be a random variable such that ~ R{k + p,Ir)- Then, according to (2.12), we 
have that for all 9 E Ir — ^, 

L,{e) = E (e<^'^^)) =e(e (e<^'^^) \k))=E (a,+,((/, - 0)"^)) . 

Using (2.5) and (2.6), we can write 
L^(0) = e(a1(,^^).(/,-^)). 

= E ((A|(i^ - e)f^-'-''^-'^ . . . (a;_i(i^ - e)f^-^^-'^ (a;(/^ - 9))-''^) . 

It follows that for all 6 G — fi. 

Setting a{9) = (log(A*_,^ (7, - 0)), log(A*_,_^ (7, - 0)), ... , log(A*_,^ (7, - 9))) , then (3.26) 
becomes 

As A; has the multivariate Poisson distribution i/, then 
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According to (3.14), we can write for all e Ir — ^, 



Therefore 



,-1 / J \ , 



'^^'^-^[A^(^)) n=^'«"--"'""'^"'^. (3.27) 



On the other hand, using (2.5), we have 



r—l / 1 \ 2 



A:._,(/.-^) = n 



1 

1 



Inserting this in (3.27), we get (3.25). □ 

Theorem 3.4 The domain of the means of the natural exponential family F = F{fj,) 
generated by the mixture fi is $7. 

Proof 

Prom (3.25), we deduce that for all e = Ir — ^, 

k,i0) = E A, ((a1,^_^^^(/, - 6)) - l) + log(A:^_^(7, - 6)). (3.28) 

i=l 

As for all i G {1, . . . , r}, the map 

(fii XI — > logA*(a;) 

is differentiable on and 

V'.ix) = iP*{x))-\ (3.29) 

then, for alH G {1, . . . , r}, we have 



(3.30) 



and for r > 2, we have 

r-l 

2 



r— 1 

(log(A;^_^(x)))' = -i Y.iP^i^)r' (3.31) 

i=l 



Differentiating (3.28) and taking into account (3.30) and (3.31), we get 



i=l 



2J ' ' " 2 

(3.32) 
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where Aq = 0. 

Let 9 e Ir — and let u be the unique element of 7^"*" such that Ir — 9 = u*~^{Ir)- Then, 
for alH G {1, . . . , r}, we have 

- ^)r^ = (i^*(u*-i(7,)))"' = (i^* (k/,))-i))"\ 

According to (2.7), this implies that for all i G {1, . . . , r}. 



r 



{p:{ir-e)r' = u\ y: • (3-33) 

\j=r-i+l J 

On the other hand, using (2.8), we can write for all i G {1, . . . ,r}, 

AUMr -0) = Ale, («*-^(^.)) = ^elWr)) = Ae._.,, (n(/.)) = 

This with (2.9) imply that for alH G {1, . . . , r}, 

A*_,^ = «2_.^i, (3.34) 

where for all i G {1, . . . ,r}, Uj arc defined in (2.10). 
Using (3.33) and (3.34), we deduce from (3.32) that 

i=l ^ ^ \j=r-i-\-\ I \i=\ 

This after a standard calculation, gives 



^^^) = E + V") ""^"^^ = ^ (E ' (3-35) 

where 

ai{Q) = '—^ + \^\. (3.36) 

As the are strictly positive, we deduce that 

fc;(e(M)) = fc;(J, - 0) C (3.37) 

Conversely, consider ?/ G O, then using the Cholesky decomposition, there exists a unique 
w G Tj"^ such that 



Ki=l / \ \i=l V IV / / \j=l 



^ 1 

Let 9 such that Ir — 9 = u*~^ (Ir) G where u = w ^ , , , &. Then 

y = u ^^ai(0)ci^ . 
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This, using (3.35), gives 
Therefore 



y = k'^{e)^k'^{Ir-^). 



n c k'^ii, -n) = k'^{Q{ix)). 

This with (3.37) imply that the domain of the means of the NEF F = F(/x) is 

Mf = k'^{Q{ii)) = n. □ 

The following theorem gives the variance function of the natural exponential family F = 
F{n). 

Theorem 3.5 For all m eft, 

>'H™) = -5p(|:j^((p;-„(™-'"" 

r 



(p;_,(m-i))~')j 



A,._,+i A,.,._,.^^(m) Xr-, A,.,._,(/n) ^ 1 



+E 



, , bj(in) 



6r_i+i(m) K-i{m) 2 

(p*{m-')y'-(PU{m-')y' 



f^i {hr-i+i{m)f 



(3.38) 



where for a// z € {1, ... , r}, hi{m) = — h \ 



Usually, for the calculation of the variance function, we set m = k'^{9) and we determine 
its reciprocal 9 = il)^{m). This is difficult to do in the present situation, however, we 
arc able to determine {P*{Ir — '4^ ■, A!_g. (/^ — ■^^(m)), and ai{ipn{m)) where Oj is 
defined in (3.36). 

Theorem 3.6 For all i e {1, . . . , r}, 



i) ai{i;^{m)) 



i-1 fi-l 



+ 



4 

r 



II) 



Ae,_,+i(m) 



+ AjAe.(m 
1 



(p;_,+i(m-i))"'-(p;_,.(m-i))" 



in) Alg.(Ir - V/.(m)) 



Proof 



(3.39) 



i) As from Theorem 3.4, m = ^^(6*) is in Q,, there exists a unique v G T^"*" such that 
m = v{Ir)- According to (3.35), we have 



v{Ir) = m = ui^ ai{tp^{m))ci j = w ( P ( ^ ^ai(V'^(m)) q I f E 



vi=l 
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Therefore 



v{Ir) =U\^P l^Yl \/^iWVi(H) Cj J (Ir) 

It follows that 



^1=1 



Or equivalently, 



1=1 



1 



u = vJ2 I (3.40) 



On the other hand, using (2.9), (3.36) becomes 

i — 1 

ai{ipn{m)) = + AjAe.(u(Ir)). 
Then using (3.40), we can write 



h AjAe- ( V — — ; -T Cj V* 

2 ' ^^ajiMm)) ' 



i-1 , . / 1 



Therefore aj('0ju(m)) satisfies the equation 



where is defined in (2.10). 

As ai{ipn{m)) > 0, we deduce that 



On the other hand, since m = v{Ir), then using (2.9), we have that 



vf = Ae,(m). 



Consequently, for alH G {1, . . . , r}. 



a,(V'^(m)) = ^ + ^ (^i^^ + XiAe^m). 
a) With the notations used above, we can write for all i G {1, . . . , r}. 
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{P*{Ir - Mm)))-' = {Pr{u*-Hlr)))-' . 



E 



ij=r-i+l 



V 



1 



1 yjai{ipij.{m)) j ^^=^-1+1 
1 

3=r-i+l 



E 9 E 



ii=l 



I M 

^aiiipi^im)) 



C; V 



Co V . 



= V 



E 



\j=i — i+l 



ajii^nirn)) 



Thus 



^ 1 



(i^*(/.-VMM)r = y „ vie,). 



(3.42) 



As m = then for all j e {1,... , r}, we have 



v{cj) = ( Ec, - y cj = (P,Vi(m-i))"' _ 



Inserting this in (3.42), we deduce that 

{P*iIr-Mm)))-'= E n\ 
Consequently 

{P*{Ir - Mm)))-'-{PU{Ir - Mm)))-' = 



(p;_,+i(m-i))"'-(p;_,.(m-^))"' 



ar-i+i{i^t^im)) 
in) According to (2.5) we have 

A*ilr - Mm)) = ^*iu*-\lr)) = A*^^iu*-\lr)) 



(3.43) 

p*{m-^)y'-(pt_,{m-^)y' 

(3.44) 



where = ^ e^. Using (2.8) and (3.40), we deduce that 

A*{Ir-Mm))=A^^^^{u{Ir)). 

= A^^*{uu*). 



A_^. UE 



} ^ I = Cj 1 1 y , = Cj 

i=i V%(V'M("i)) / Vi=iv%(^MM) 



c,- V 
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= A_ij* I V — — Cj V* 



Therefore 



Thus, using (2.4), we obtain 



A*{Ir-Mm)) 



It follows that for alH G {1, . . . , r}. 



1 



9 



1 



Al{Ir-Mm)) 
Using (2.9), we deduce that for all i G {1, . . . ,r}. 



^'-^l"lga,(V^,\m))^^ 



AUMr-^rn)) 



ar-i+i(V'//(m))' 
Therefore 



A_^.{Ir - tpij.{m)) - 



ar-i+i{ipiu{m))Ar-i{m) ' 

^e,._i+i(m) 



ar-i+i{tp,j,{m))' 



□ 



We are now in position to give the variance function of the natural exponential family F 
stated in Theorem 3.5. 

Proof of Theorem 3.5 

We have that for all m G Mp = ^, 

Vpim) = k^iMm)). 

Differentiating (3.32) and using (3.30) and the fact that for alH G {1, . . . , r} and x G O, 

{{P*{x))-')' = -P{{P*{x))-'), 

we get for all 9 & Ir — ^, 
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k';,{e) = --p({p:{ir-d))-') 

+ (A.-i+lAl,^(/, -6)- Xr-^A*_,^^^{Ir - 0) + \) {P {{Pt i^r - 0))-')) 
i=l 

- ^ Xr-^AU,^,iIr - 6) {{P^lilr - 0))-' - {P*{Ir - 9))-') ® {P* {Ir - 9))-' . 
i=l 

It follows that for all 9 e Ir — ^, 

k';,{9) = -\p({p;{ir-9))-') 

i=l 

X [{iP*{I, - 9))-' - {PUilr - e))-') ® {{P*{Ir - e))-' - {PUilr ' 0))-')' ■ 

We need only to replace 9 by ip^{m), then insert (3.39), (3.43) and (3.44) to get the 
expression of the variance function of F = F{iJ,) given in (3.38). □ 
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